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Abstract 

We review and discuss the equations governing the distribution of work done on a system which is driven out of 
equilibrium by external manipulation, as well as those governing the entropy flow to a reservoir in a nonequilibrium 
system. We take advantage of these equations to investigate the path phase transition in a manipulated mean- 
field Ising model and the large-deviation function for the heat flow in the asymmetric exclusion process with 
periodically varying transition probabilities. To cite this article: A. Imparato, L. Peliti, C. R. Physique 6 (2005). 

Resume 

Distributions du travail et de la chaleur dans des systemes hors equilibre. Nous passons en revue et 
discutons les equations regissant la distribution du travail effectue sur un systeme manipule hors d'equilibre, ainsi 
que celles qui regissent le flux d'entropie vers un reservoir dans un systeme hors d'equilibre. Nous exploitons ces 
equations dans l'etude de la transition de phase dans les chemins d'un modele d'Ising champ moyen manipule et 
de la fonction des grandes deviations pour le flux d'entropie dans le modele d'exclusion asymetrique a probabilites 
de transition periodiques dans le temps. Pour citer cet article : A. Imparato, L. Peliti, C. R. Physique 6 (2005). 
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Mots-cles : Proces hors equilibre; Distribution du travail; Flux d'entropie 



1. Introduction 



In the recent years we have seen an important outburst of activity in the field of nonequilibrium thermody- 
namics and statistical mechanics, sparked by the discovery of a number of results of remarkable generality and 
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impact [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18]. Some of these results yield predictions on the properties of 
the distribution of the work performed on a system as it is manipulated, while others describe properties of the 
distribution of the entropy generated in the nonequilibrium process. It is clear that in order to make good use of 
this information, it is advisable to investigate in detail the properties of the distribution of both work and heat 
in a manipulated system. 

In the present Contribution we discuss the work and the heat distribution in some systems, pointing out some 
aspects that we have found interesting enough to be brought to the attention of our colleagues. We first consider 
the expression of the generating function for the work distribution in a "large" system. Although our results 
formally hold in the limit of infinite size, they are of interest also for systems small enough that the relevant 
energy barriers are of the order of a few ksT, well within the range in which work fluctuations can be observed. 
In a large system, as already pointed out by Ritort [19] and discussed in [20], the generating function of the work 
distribution is dominated by the contribution of phase space paths which satisfy an ordinary differential equation 
akin to a classical equation of motion. We have found that in some situations these paths can exhibit a singularity, 
for some protocols, as a function of the variable conjugate to the accumulated work. This singularity is similar to 
a phase transition taking place in the space of paths. We shall exhibit a simple model in which this phenomenon 
takes place and attempt to define the corresponding phase diagram. 

We also consider the heat flow distribution in a general Markov process, whose differential equations were 
derived by Lebowitz and Spohn [8]. We were able to evaluate the solution of this equation for the asymmetric 
exclusion model with periodically varying transition probabilities, and to exhibit that the Gallavotti-Cohen [5] 
symmetry also holds in this case. 



2. The generating function for the work distribution 



We first briefly review the derivation of the generating function for the work distribution of manipulated systems, 
which may be found in [20]. Let us consider a system with discrete states i of energy Hi(n), where \jl is a parameter 
which is manipulated according to some protocol fJ,(t), starting at t = 0. The evolution of the system is described 
by a markovian stochastic process: given, for any two states i, j, the transition rate kij(t) from state j to state i 
at time t, the system satisfies the set of differential equations 

E [M*)Pi(*)-M*)Pi(*)]. (i) 

where Pi(t) is the probability that the system is found at state i at time t. Let p° q (A0 represent the canonical 
equilibrium distribution corresponding to a given value of fi. We have 

p -l3H t { M ) 

P?M = , (2) 



where = J2 i e t3H *^ = e 0F » is the partition function corresponding to the value /i of the parameter, and F t 



the relative free energy. The transition rates fc»j(t) are compatible with the equilibrium distribution p^ q (/i), i.e., 
one has, for all i, 

E [M*)pr(M*)) - M*K q (M*))] = o. (3) 

We assume that the system is at equilibrium at t = 0, and therefore, that Pi(t) satisfies the initial condition 

P^=0)=P- q (M(0)). (4) 

We now consider the joint probability distribution $i(W, t) that the system is found in state i, having received a 
work W, at time t. The set of differential equations satisfied by the distribution functions $i(W, t) reads 

= E \kiAmAw,t) - M*)*«Wt)] - (5) 
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The joint probability distribution $i{W,i) satisfies the initial condition *i(W,0) = 8{W) p° q (/i(0)). Then the 
state-independent work probability distribution P(W,t) is defined by P(W,t) — £\ t)- It is convenient to 

introduce the generating function of $i with respect to the work distribution, defined by 

#i(A,t) = J AW e xw ^i{W,t). (6) 

We assume that $j(W, f) vanishes fast enough, as \W\ — > oo, for >]/i(A,i) to exist for any A. The function V&i 
satisfies, for all A, the initial condition \I/i(A, to) = p" q (/i(0)), and evolves according to the differential equation 

ft* i (A,t)= ]T [fcij^j - fcji^i] + A/z ^(A,t). (7) 

Exploiting (3), it is easy to verify that, if A = —f3, the solution of (7) with its initial condition reads 

**(-/?, *) = —7 — = p- q (MW). (8) 

We can thus straightforwardly check that the solution of (7) verifies the Jarzynski equality [6]: 

e -0W\ = Z »{t) = e -/3(F( M (t))-F( M (0)))_ , q n 
/ Z M (0) 

It is thus possible, in principle, to evaluate the probability distribution function of the work W by solving the 
equations (5) or (7) for all the microscopic states i. This approach has been implemented in [21] for a simple 
model of a biopolymer. 

The approach discussed in the previous section becomes quickly unwieldy as the complexity of the system 
increases: the dimension of the system (5) is equal to the number of microscopic states of the system. When the 
system considered is characterized by a large number of degrees of freedom, it is convenient to introduce some 
collective variables, and an effective free energy, in order to reduce the complexity of the problem. The assumption 
underlying this approach is that the system reaches on a comparatively short time scale a quasicquilibrium state 
constrained by the instantaneous value of the collective coordinate. Thus, on the the time scale of the experiment, 
the state of the system can be well summarized by the collective coordinate, with the corresponding free energy 
playing the role of the hamiltonian. 

Let us now consider a system characterized by a generic equilibrium free energy function J-^{M), where \i 
is again the parameter which is manipulated, and M is some collective (mean-field) variable. We assume that 
the system dynamics is stochastic and markovian: let P(M, t) denote the probability distribution function of the 
variable M at time t, then its time evolution is described by the differential equation 

dP 

^ = np, (io) 

where "H is a differential operator which depends on the parameter fi. We require that the operator Ti is compatible 
with the equilibrium distribution function of the system, i.e., that the relation 

He-^ (M) =0 (11) 

holds for any value of ju. 

The work done on a system during the manipulation, along a given stochastic trajectory M(t), is given by 



t 

W = 

o 



(12) 



Using the same arguments as for the discrete case, one finds that the time evolution of the joint probability 
distribution <I>(M, W, t) of M and W is described by the differential equation 

It can be easily shown that the solution of (13) identically satisfies the Jarzynski equality [22]. 
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Equation (13) becomes much easier to treat if one introduces the generating function <f(M, A,t) for the work 
distribution: 

#(M,A,t) = jdWe xw ${M, W,t). (14) 
Equation (13) thus becomes 

with the initial condition *(M, A, 0) = p^ 0) (M) = e-^C) (M) /^(o)- 

As shown in [22], one can derive a path integral representation of the solution of (15), taking for the differential 
operator 7i the expression 

oo _ fc 

=Eiw«)i as) 

fc=0 

(The coefficients gk(M) also depend on fi, but this dependence is understood to lighten the notation.) One obtains 
the formal solution 

M(t f )=M 

#(M,A,t f ) = J dM J V-fVM exp{S[ 7 ,M]} *(M ,A,0), (17) 

M(0)=M 

where S[7, M] = J Q ! dt C(t). The "lagrangian" £ is given by 

C(t)= Lm + H( 7 ,M) + \i1^) , (18) 

V °f I / 7 (t),M(t), M (t) 

where the "hamiltonian" 7i(7, M) is defined by 

W( 7 ,Af)=^ 7 fc 5fe (M). (19) 

fc=0 

Let iV indicate the size of the system, and let us define the "intensive quantity" m = M/N . We can thus define, 
in the thermodynamic limit N — > oo, m — const., the densities 

Mm) = Jim ff(7 , m) = ffin /(t) = lim M (20) 

iV^oo iv iv^oo iv N^oo iV 

with ^(t) = 7m + i?( 7 , m) + Xjldf^/dfj,. As discussed in [19,22], when the system size iV is large enough, the path 
integral in (17) is dominated by the classical path (7c (t), m c (t)), solution of the equations 

SS dH SS . flff,.. a 2 /, 



5l {t)-^ rh = -^' ^ =0 ^ = ^ + AA ^? (21) 

We shall now see that by requiring the system is in equilibrium before the manipulation starts, we impose an initial 
condition on these equations. In order to evaluate the integral over Mo in (17) with the saddle-point method, we 
note that ~$>(M, A, 0) appearing on its rhs, is given by p^L,(M). Furthermore, from the definition of £(t), (20), it 
follows that 



j dt£(t) = m tf j tf - m 7o + J dt [—ym + H + A/i<9 M / M ] . 



(22) 



Thus, substituting (22) into (17), and taking the derivative with respect to mo = Mo/N, we obtain the saddle-point 
condition 



(23) 



In this way one can devise a strategy to evaluate ^(M, A,if) for a given manipulation protocol fi(t), when the 
system size N is large enough. One has to solve the classical evolution equations (21) with a two-point boundary 
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condition: namely, (23) should be imposed at t = 0, and the condition Nm(tf) = M should be imposed at 
the final time t{. Once the relevant classical path (-y c (t),m c (t)) has been evaluated, one can obtain the action 
density s[~/ c ,m c ] = limjv^oo <S[7c, Nm c ]/N from the expression S [ 7c ,m c ] = f Q l dt e(t). We are interested in the 
state-independent work probability distribution 

P(W,t { ) = J d\e- xw T(\,t f ), (24) 

where we have defined r(A,tf) = J AM ^>(M, X,t f ). We shall now see that evaluating r(A,t f ) identifies a well- 
defined boundary condition on y c (tf)- We have indeed 

M(t f )=M 

#(M ,A,0). (25) 



r(A,t f ) = J dMdMo J V-/VM cxp N J Atl(t) 



M(0)=M 

In order to evaluate the integral over M with the saddle point method, we notice that, upon derivation of the rhs 
of (22) with respect to m tf , we obtain the condition 7f = j(t{) = 0. Thus, the equation of motions (21) have to 
be solved with the initial and the final conditions that we have derived: let (7* (£), m*(t)) denote the solution of 
equations (21) satisfying these conditions. For each value of A, taking into account its initial condition, we obtain 
the following saddle point estimation for T(X,tf): 

r(A,* f )oc eXpl ^ (A)] , (26) 

where 

tf 

g{\) = J dtll^-pf^ml). (27) 


In this equation, £*(t) is £(t) evaluated along the classical path (7*(t), m*(t)). In order to evaluate the integral on 
the rhs of (24), we use the saddle point method again, and obtain 

P{Nw,t t ) = Nexp{N[-X*(w)w + g(X*(w))]}, with g'{X*)=w, (28) 

where M is a normalization constant. Notice that the saddle point estimate for P(W,t{) obtained in this way, 
implies that the distribution becomes more and more sharply peaked around its maximum value as N — > 00. This 
is compatible with the expectation that the work fluctuations becomes relatively smaller as the size of the system 
increases, and in the limit iV — > 00, which can be thought as the limit of a macroscopic system, no work fluctuations 
are observed, and the work done on the system during the manipulation takes one single value, corresponding to 
the most probable value of P(W,t{). 



3. A path phase transition 



In this section we consider an Ising model in mean-field approximation, which evolves according to a Fokkcr- 
Planck equation, whose differential operator reads 



H- =uj N 



d 



dM 

leading to the hamiltonian 



( — ) 

\dM J 



+ : 



dM 



H (7, m) — u>o 



\dm ) + 



0-1 2 

13 7 



where the free energy density f(m) = T(Nm)/N is given by 



f(m) = - '-^n - hm + /3 1 



1 + m 



log 



1 + m 



+ 



1 — m 



log 



1 — m 



(29) 



(30) 



(31) 
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We take the magnetic field h(t) as the external parameter that varies with the time and drives the system out 
of equilibrium. The equations of motion (21) thus become 



dH df n , „ 

m = = -wot: 2fc B Two7; 

cry dm 



dH _ d 2 f d 2 f ■ 

a ^ V a a = Wotj— 9 7 - Aft, 



(32) 



The magnetic field /i(£) is taken to vary according to the linear protocol 



h(t) = h + {hi - h )- 



0<t<t { . 



(33) 



Here we consider the case of Ising model below the transition temperature, i.e., J — 1.1, and the initial and final 
value of the magnetic field h(t) are taken to be ho = —hi = — 1. In the present section we set (3 = 1/ksT = 1. We 
plot in figure 1 the probability distribution of the work done on the model, as obtained from (28), for two values of 
the manipulation rate. In the same figure the histograms of the work obtained by simulating the process are plotted. 
The process is simulated by integrating the corresponding Langevin equation, using the Heun algorithm [23,24]. 
The agreement with the curves as obtained from (28) is satisfactory. In the insets of the same figure, we plot the 
quantity P(w) defined by 



P(w) = exp [-/3Nw] P(w). 



(34) 



On the one hand we find J dwP(w) = exp [—/3AF] = 1 as predicted by the JE (9), while on the other hand the 
histogram obtained by the simulations exhibits no point (no realization of the process) with w < = AF. Thus 
the work distribution obtained by the simulation of the process cannot reliably be used for estimating AF. This 
is a typical example of how the lack of knowledge of the tails of the work distributions in micro-manipulations 
experiments hinders the possibility of using (9) to evaluate free energy differences. 




a, 




Figure 1. Probability distribution function P(w) for the system described by the differential operator (29) with equilibrium free 
energy (31), manipulated according to the protocol (33), with J\ — 1.1, ho — —hi — — 1, and tf — 2 (left panel), tf = 4 (right 
panel). Continuous line: probability density P(w) of the work "per spin" w — W/N, with N — 100. The histogram of the work is 
obtained by 10000 simulations of the process, see text. Insets: P(w) as given by (34), whose integral verifies the Jarzynski equality. 

We now consider the classical paths m*(t, A), solutions of equations (32), and thus contributing to the probability 
distribution P(w) via (28). In the left panel of figure 2, we plot m* as a function of t for different values of A, 
obtained by numerical solution of equations (32), for U = 2: we observe that the trajectory m*(i,A) varies 
continuosly as A is varied. As a consequence, since the work done on the system along each trajectory m*(t,A) 
reads w — — dt'h(t')m*(t' , A), w turns out to be a continuous function of A (figure 2, right panel). Furthermore, 
since w and the saddle point value A* are related by (28), the function g(X) as given by (27) is a differentiable 
function with respect to A, as shown in the inset of figure 2, right panel. We now consider a faster protocol, tf = 0.2: 
the results are plotted in figure 3. One can clearly see that the classical paths m*(t, A) exhibit a discontinuity for 
A = 0.5, jumping from negative to positive values. Accordingly, w(X*) exhibits a discontinuity at A* = 0.5, as 
shown in the right panel of figure 3. This is reflected in the appearence of a cusp in the function p(A), at A = 1/2, 
see the inset of figure 3. We find that for any value of r = (hi — ho)/U it is always possible to find a value of 
J = J*(r) such that for J > J*(r) the function w(X*) cxibits a discontinuity at A = 1/2, and thus the trajectories 
m*(t, A) exhibit a path phase separation. In figure 4, left panel, we show the system phase diagram: in the upper 
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0.5 1 1.5 2 -3 -2 -1 1 2 3 4 

t X* 

Figure 2. Left: plot of m* as a function of t for different values of A, with J — 1.1, ho — —hi — — 1, and if — 2. The values of A 
vary between A — —5 (bottom curve) and A — 5 (top curve), with a step AA — 0.2. Right: plot of w as a function of A* , as defined 
by (28). Inset: plot of g as a function of A as given by (27). 



1 




-0.5 



-1 1 ' *— ■ 
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Figure 3. Left: Plot of m* as a function of t for different values 
A vary between A — —5 (bottom curve) and A — 5 (top curve), 
function of A*, as defined by (28). Inset: plot of g as a function < 
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f A, with J = 1.1, h = -hx = -1, and i f = 0.2. The values of 
ith a step AA — 0.2. Thick line: A — 0.5. Right: plot of w as a 
A as defined by (27). 




Figure 4. Left: Ising model phase diagram, in the r — (hi — ho)/t{, J plane. In the upper right part of the diagram the system 
exhibits path phase separation, while in the lower left part no phase separation is found. The points correspond to the value J* (r) 
such that for J > J* (r) the function u>(A*) cxibits a discontinuity at A — 1/2. Right: Path Helmholtz free energy for the system 
described by the differential operator (29), manipulated according to the protocol (33). Full line: J — 1.1, h — —hi — —1, and 
t{ — 0.2. Dotted line: J — 1.1, ho — —hi — —1, and t{ — 2. The curves are arbitrarily shifted for clarity. 



right part of the diagram the system exhibits path phase separation, while in the lower left part the trajectories 
m* (£, A) vary continuously as A is varied. 
Let us now define the function 
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4>{w) = - lim —logP(Nw,t c ). (35) 

N — >oo iv 

The functions g(X) and <j>(w) are related by a Legendre transformation: 

4>(w) = inf (-p(A) + Aw) = -g(A*(w)) + A*(w) w, (36) 

where A*(w) is the solution of (28). Thus A and w act like thermodynamically conjugate variables. These functions 
can be thus interpreted in terms of path thermodynamics: <?(A) can be viewed as a path Gibbs free energy, while 
4>{w) is the corresponding Helmholtz free energy. This analogy was first pointed out in [19,20]. 

How do we interpret the singularity in g(X)? In strict analogy with thermodynamics, the discontinuity of g'(X) 
with respect to its independent variable A corresponds to the linear behavior of its Legendre transform (j>(w) 
between (w+,4>{w+)) and (w-,cf>(w-)), where w± are the values of w either side of the discontinuity. Thus the 
"path phase coexistence" appears as an exponential dependence of P(Nw,tf) = exp(— N<f>(w)) on W in a certain 
interval. The presence of such exponential tails was conjectured in [19] for a system of independent spins, but we 
were only able to exhibit them in an interacting system like the present one. 



4. The distribution of heat flow in a Markov process 



In this section we discuss the equation governing the time evolution of the probability distribution function of 
the entropy which flows into the enviroment surrounding a stochastic system which evolves across its phase space. 
We assume that the system at issue has a discrete phase space and its time evolution is a stochastic markovian 
process described by (1). For simplicity, we consider a stochastic dynamics with a discrete small time scale r, 
such that the jumps between states take place at discrete times tk = kr. We consider a generic path lo defined by 
ui(t) — ik iff tk < t < tk+i, with k — 0,1, ... , M, with tui+i = U, and define the time-reversed path lo by uj(t) = ik 
for tk+i <t<t k , where t = to + U — t. Let us define the quantity Q(u>) by 



QH = - In 



E ln 



(37) 



where V(lo) is the probability of the forward path lo (conditioned by its initial state io) and V(u>) is the probability 
of the time-reversed path lo, conditioned by its initial state im = if and subject to the time-reversed protocol 
Kij(t) = Kij(t) [8,9,10,16]. We have assumed that, if K XJ {t) > at any time t, one also has K^t) > 0. 

It is worth noting that, if the detailed balance conditions holds for the transition rates Kij(t), and the energy 
Hi{t) is associated to the state i of the system, we have Kji(t) /Kij(t) = oxp{[Hi(t) — Hj(t)] /T}, and thus 
Tin [Kji(t)/Kij(t)] represents the heat exchanged with the reservoir in the jump from state j to state i. (In this 
section we set fee = 1.) Thus the quantity Q(ui), defined by (37), is the entropy which flows into the reservoir as 
the system evolves along the path w [8,14,15,25]. Let us define As^ as the entropy which flows into the reservoir 
as a result of the jump of the system from state j to state i Astj = log [Kji(t)/ Kij(t)]. The differential equation 
governing the time evolution of the joint probability distribution function (<?,£) reads [25] 



d$i(Q,t) 
dt 



E 



(-A Sij ) n d n ^(Q,t) 



Kji^i(Q,t) 



(38) 



By introducing, for each i, the generating function ^>i(X,t) = J dQ exp(AQ) $i(Q,t), and taking into account the 
expression of Asij, we obtain the master equation 



dt 



E 



= ^^(A)* 3 (A,t). 



(39) 



which was first derived by Lebowitz and Spohn in [8]. In the case of time-independent transition rates Kij, or 
of transition rates which depend periodically on the time, it can be useful, in order to evaluate the distribution 
function $(Q,t) = 5^^<&(Q,t), to introduce the large-deviation function. In the long-time limit, the generating 
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function ^(A,£) = i) is dominated by the maximum eigenvalue g(X) of the matrix H(A) = (Hij(X)), 

which appears in the master equation (39). Therefore, we have, for long times t, 

®{\,t)<xexp[tg(X)]. (40) 

By using the last equation and using the definition of the generating function, one obtains the probability distri- 
bution of the entropy flow in the long time limit: 

*(Q >t ) = f ^ c - A «*(A,t) oc e t9 ( A *)- A * Q , (41) 

J 27T1 

where A* is the saddle point value implicitly defined by dg/d\\ x , = Q/t. If we introduce the entropy flow per 
unit time q = Q/t, we obtain the large-deviation function 

f(q) = g(X*) -X*q = lim ) log <t>(tq, t). (42) 

Note that the functions g(X) and f(q) are Legendre transform of each other, and can be then interpreted in 
terms of path thermodynamics: g (A) can be viewed as a path Gibbs free energy, while f(q) is the corresponding 
Helmholtz free energy. The connection between the generating function \&(A) and the thermodynamic formalism 
for dynamical systems [26] has been investigated in Ref. [27]. 

If the system is characterized by a small number of states, one can explicitly solve the equations (39) , and thus 
obtain the total generating function ^(A,i) ee - "fy i (X,t) — (e A< ^) . While on the one hand this direct approach 
becomes rapidly impracticable, as the system phase space size increases, on the other hand (39) suggests a practical 
computational scheme to evaluate the generating function \P(A,t). Since ^(X,t) = f Vut V(u>t) e A<3 ^*- ) , we have 

^l*)=(0) A *(A,t), (43) 

where {■ ■ ■} x is the average in the weighted ensemble V(iOt) exp [XQ(u>t)] /^{X, t), where ^(A,^ = Z\ is the 
"partition function" of this ensemble, which will be called the "A-ensemble" in the following. The solution of (43) 
thus reads 



#(A,f) = exp 



dA' (Q) y 



(44) 



Following [25], in the present paper we consider a procedure to evaluate Vl'(A,t) which generates trajectories in 
a suitable entropy-flow weighted ensemble. The direct simulation of trajectories in the A-ensemble is hindered by 
the fact that one should already know the exact expression of the ensemble partition function, i.e., the function 
*(A,t), which is the unknown quantity at issue. 

To avoid the problem of the direct evaluation of *(A, t), following [28], we introduce a generic functional of the 
paths n(w), and write 



J Vujt (QM/n(^))n(^)?(^)e A «"'» = (Q/n) A[n 
fVuj t (i/n(u t ))ii{Lo t )v(Lo t )e^^) - <i/n) An 



(Q/X = r-r... Ti TttZTTvttZ . v^7-7T73^wi = TTTm > ( 45 ) 



where (. . .) A n indicates the average in the new P(wt)n(wt) exp [AQ(wt)] ensemble, which will be indicated as the 
(A, II)-ensemble in the following. 

We choose the functional of the path II (u) as discussed in [25]. The probability of a given path lj reads 

V{ljj) = Cli N ,i N _ 1 {tN-l)Qi N _ 1 ,i N _2 

(tjv- 2 )...fii 1 ,i (to)Pio> (46) 
where the transition probabilities £li,j(t) are defined as Qij = rK ij (t), and = 1 — ^2^^ Slji(t). We now define 
the new transition probabilities Q,i,j = rKij (Kji/Kij) x , and f2j,j = 1 — ^i*> an d choose the functional 

II(w), such that [25] 

M , .. . . . 

nH = n n — (^), with ^(t) = U j{t)/ n jj{ t), Vit]: ( 47 ) 

k — 1 
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Recalling the definition of Q(ui), (37), we obtain that the probability in the (A, il)-ensemble is given by 

V(u)U(u) exp [AQH] = n tN , iN _ 1 ?l iN4N _ 1 . . .Q. ilAo p° , (48) 

and thus {Q) x can be evaluated by using (45), for the particular choice of II, as given by (47). Note that (48) 
implies that the one can generate a trajectory in the (A, n)-ensemble by simply simulating the process with the 
Qij(t) transition probabilities. 

In [25], the the feasibility of the method was illustred by applying it to a nonequilibrium system characterized 
by a large phase space, and evolving according to a stochastic dynamics, namely the simple asymmetric exclusion 
process (ASEP) [29]. Such a system consists in a one-dimensional lattice gas on a lattice of L sites. Each site of the 
model is either empty or occupied by at most one particle. Each particle can jump into an empty nearest neighbor 
site with transition rates per time unit K + (rightward) and K- (leftward). The system is kept in an out-of- 
equilibrium steady state since its first and last site are in contact with two particle reservoirs, at densities pA and 
pB respectively. By taking pA> Pb and K+ > K-, one observes a net particle current from the left to the right 
reservoir. In [25] we considered an ASEP model with constant parameters pA, Pb, K+, K~. Here we consider a 
system with time-dependent parameters. We take K + — 1, K- — K°_ (1 + sin(27rf)) + e, with K°_ = 0.75, e = 10~ 3 , 
a = 1, 7 = 0.27, (3 = 1, 5 = 0.27, where a-y are the rates of jump from-into the left A reservoir, respectively, 
and /3-8 are the rates of jump into-from the right B reservoir, respectively. With this choice of parameters, at 
t — 2kir, we have pA = 0.75, pb = 0.25 [30]. The duration of a single trajectory is taken to be U = 10 • 27r. A 




A q 

Figure 5. Left: Plot of g(\) as obtained by combining (40) and (44), for the ASEP model. The function g(X) vanishes for A — 0. 1 
which corresponds to the normalization condition and to Seifert's fluctuation relation [16] respectively. Right: Histogram of the 
entropy flow per time unit q, corresponding to 1000 unbiased trajectories. Full line: probability distribution function <&(tfg, tf), 
obtained with the trajectory sampling algorithm. 



direct evaluation of the function g{\) by solving the 2 100 equations (39) is out of question. We thus apply our 
trajectory simulation approach to the ASEP model. We consider trajectories with elementary time step r = 0.01: 
at each time the transition probability between two states is given by the transition matrix Kij(t). For each value 
of A we generate Af = 10000 sample trajectories and calculate the entropy flow Q, as defined by (37), for each 
trajectory. Then, for the given value of A, by averaging over the Af trajectories, we compute the quantity {Q) x 
using (45). Finally, by combining (40) and (44), we obtain the function <?(A) which governs the long time behavior 
of *(A, t). This function is plotted in figure 5, left panel. It can be seen that it vanishes at A = 0, 1 and is symmetric 
with respect to A = 1/2. The fact that g(0) = corresponds trivially to the normalization condition over all the 
possible trajectories. On the other hand, the fact that the function g vanishes at A = 1, is a non trivial result, and 
corresponds to Seifert's fluctuation theorem [16]. The symmetry around A = 1/2 corresponds to the Gallavotti- 
Cohen fluctuation relation. We are now able to calculate the large-deviation function /(g) defined in (42). We 
check as follows that the quantity f(q) actually gives the entropy distribution function $(g, t) oc exp [t /(g)] for 
the present model. We simulate the unbiased diffusion process by using the transition matrix Kij(t), and measure 
the entropy flow along 1000 trajectories. We then plot the histogram of the measured entropy flow per time unit, 
together with the function exp [t /(g)], see figure 5, right panel. The agreement between the histogram and the 
predicted entropy distribution $(g, t) is excellent. 
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5. Discussion 



We have seen that the equations governing the evolution of the work and heat flow distributions in out-of- 
equilibrium systems can be harnassed to yield interesting information. We have analyzed the generating function 
of the work distribution in "large" manipulated systems and shown that it may exhibit a "path phase transition" , 
which corresponds to the presence of an exponential behavior in some work interval. On the other hand, we have 
shown that it is possible to evaluate the large deviation function for the entropy flow in a stochastic process via 
a biased simulation technique, and applied the method to exhibit the Gallavotti-Cohen symmetry in an out-of- 
equilibrium system with periodically varying parameters. A similar result has been recently obtained by Ge and 
collaborators [31]. 
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